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Abstract. We have determined the mass function of loose groups of galaxies in the Las Campanas Redshift Survey. 
Loose groups of galaxies in the LCRS range in mass from M ~ 10 12 Mq to 10 15 Mq. We find that the sample is 
almost complete for masses in the interval 5 • 10 13 — 8 • 10 14 Mq . Comparison of the observed mass function with 
theoretical mass functions obtained from N-body simulations shows good agreement with a CDM model with the 
parameters fi m = 0.3, Qa = 0.7 and the amplitude of perturbations about as = 0.78 — 0.87. For smaller masses 
the mass function of LCRS loose groups flattens out, differing considerably from the group mass function found 
by Girardi and Giuricin (2000) and from mass functions obtained by numerical simulations. 
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1. Introduction 

One important constraint of cosmological models is pro- 
vided by the mass function of clusters and groups of galax- 
ies. In the case of popular ACDM models the amplitude 
and shape of the mass function depend primarily on the 
mean density of matter in the Universe, fl m = O c + Oj, , 
where f2 c and Qb are the mean densities of the cold dark 
matter and baryonic matter in units of the critical cos- 
mological density, respectively. The amplitude of the mass 
function also depends on the amplitude of the power spec- 
trum of density fluctuations, which can be characterised 
by the as parameter (the linearly extrapolated rms den- 
sity fluctuations in a sphere of of 8 h~ x Mpc radius). The 
local abundance of rich clusters of galax ies can be used 
to estimate o% (Bah call, F an & Cen 1997 , de Theije, van 
Kampen & Slijkhuis [l998] Cen |l998j ). 

Determination of the cluster mass function consists of 
two tasks: calculation of cluster masses, and estimation of 
the spatial density of clusters. Masses of clusters can be 
derived using either X-ray data and the mass-temperature 
relation or data from optical surveys using the velocity 
dispersion of galaxies in clusters ( virial masses) . In a pio- 
neering study by Bahcall and Cen ( 1993| ; hereafter BC93), 
both mass determination methods were used. Masses were 
derived for Abell clusters of all rich ness c lasses and for 
groups of galaxies. Biviano et al. ( |1993| ) and Girardi 



et al. (1998, hereafter G98) calculated virial masses of 



flux-limited sample of galaxy clusters (HIFLUGCS, the 
Highest X-ray FLUx Galaxy Cluster Sample) to obtain a 
mass function. The spatial density of massive clusters of 
galaxies, according to G98, exceeds the density found by 
BC93 by a factor of almost ten. As shown by G98, this 
difference is mainly due to the fact that BC93 assumed a 
one-to-one correspondence between the richness of a clus- 
ter and its mass, whereas in reality this relationship has 
a large scatter. A new determination of the mass function 
of loose groups of galaxies has been provided by Girardi 



and Giuricin (2000) 



nearby clusters. Reiprich & Bohringer ( |2002| ) used a X-ray 
Send offprint requests to: P. Heinamaki 



The sample of loose groups of galaxies identified in the 
Las C ampanas Redshift Survey (LCRS) by Tucker et al. 
( 2000l hereafter TUC) provides a possibility to derive a 
new independent estimate of the mass function of clus- 
ters and groups of galaxies. Loose groups represent the 
most numerous component of galaxy clustering, and so 
we should be able to determine their statistical properties 
much better than those for rare rich clusters of galaxies. 
The virial masses of loose groups of galaxies in the LCRS 
range from M ~ 10 12 Mq to 10 15 Mq; these masses were 
estimated by TUC using velocity dispersions of galaxies 
and harmonic radii of groups. 

In this paper we shall estimate the mass function of 
LCRS loose groups. Our main task is to derive the spa- 
tial density of groups of various masses. The volume of 
the LCRS sample is well determined, thus we hope to get 
an unbiased estimate of the spatial density of groups. We 
shall compare the empirical mass function with theoret- 
ical mass functions found using numerical simulations of 
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structure evolution. In our simulations we assign to a halo 
all the particles identified as members of the halo by halo 
finder. In observations this means that the velocity dis- 
persion of galaxies is assumed to be identical to that of 
the dark matter. This comparison allows us to check the 
consistency of structure evolution models with empirical 
data and to find the set of cosmological parameters which 
brings models into agreement with data. We also deter- 
mine the mass interval where the sample of groups is not 
influenced by selection effects. 

2. The Data 

2.1. Observations 



The LCRS (Shcctman et al. 1996 ) is an optically selected 
galaxy redshift survey where a multi-object spectrograph 
was used to measure simultaneously redshifts of 50 or 112 
galaxies. Extending to a redshift of z ss 0.2, the catalogue 
covers 6 slices of size on average 1.5° x 80°, containing 
23,967 galaxies with measured redshifts within the official 
survey photometric and geometric limits. Three slices are 
located in the northern Galactic hemisphere, centred at 
declinations S = —3°, — 6°, — 12°, and the other three 
slices are located in the southern Galactic hemisphere, 
centred at declinations S = -39°, - 42°, - 45°. TUC 
applied a friends-of-friends (FoF) percolation algorithm 
to extract the catalogue of Las Campanas loose groups 
of galaxies (hereafter LCLG). The linking parameter was 
chosen to get a density enhancement limit of Sn/n > 80. 
The minimum group membership was chosen to be three. 

Because the spectroscopy was carried out for each field 
either via a 50 or a 112 fibre multiobject spectrograph, the 
selection criteria varied from field to field. The nominal 
apparent magnitude limits for 50 fibre fields were 16.0 < 
R < 17.3 and for 112 fibre fields the limits were 15.0 < 
B < 17.7. According to TUC great effort was put forth 
in accounting for these field-to-field sampling variations. 
The general properties of the 50 fibre and the 112 fibre 
groups agree well with properties of groups found from 
other surveys, thus we shall use the whole TUC group 
sample. 

Only groups with redshifts 10,000 < cz < 45,000 
km s _1 were included into the LCLG sample. TUC con- 
cluded that the LCLG is a useful sample for a variety of 
studies requiring an unbiased collection of loose groups. 
It is based on the LCRS galaxy sample, which is the first 
redshift survey that can claim to enclose a reasonably fair 
sample of the nearby universe. 

The complete LCLG list includes 1495 groups. TUC 
also introduced a "clean sample" , where groups with four 
potential bias factors are excluded: 1) groups which are 
too close to a slice edge, 2) groups which have the cross- 
ing time greater than the Hubble time, 3) groups with 
the corrected velocity dispersion zero what can happen 
since TUC subtracted an redshift measurement error of 
67 km/s in quadrature from each group velocity disper- 
sion, and 4) groups containing a 55 arc-sec orphan galaxy, 



i.e. a galaxy with no measured redshift. The last bias was 
caused by technical reasons (the fibre separation limit, 
which prevents the observation of neighbouring galaxies 
within 55 arc-seconds of each other). In the full sam- 
ple this effect was compensated for by reintroducing lost 
galaxies, assigning to them a redshift equal to that of its 
nearest neighbour, convolved with a Gaussian of width 
cr = 200 km s _1 . 



2.2. Simulations 

A low density CDM universe with a cosmological con- 
stant (ACDM) is widely regarded as the best model com- 
patible with most of the currently available data; e.g. 
with the microwave background anisotropy measured by 
BOOMERANG (de B ernardis et al. |2000D and MAXIMA 
I (Hanany et al. 2000 ) experiments, and with data on the 
large scale structure of the universe: two-point correlation 
functions and power spectra, and the spatial density of 
mass-limited samples of galaxy clusters (Goyernato et al. 



1999| , Colberg et al. gOOO] Pierpaoli et al. gOOlj) 



For the present study we employ a flat cosmologi- 
cal model (f2 TO + Q\ = 1) with the following parame- 
ters: the matter density f2 m = 0.3, the baryonic density 
fib = 0.04, the vacuum energy density (cosmological con- 
stant) Sl\ — 0.7, and the Hubble constant h = 0.7 (here 
and throughout this paper h is the present-day Hubble 
constant in units of 100 km s _1 Mpc -1 ). The simula- 



tions were performed using a P3M code (Couchman 1991 ) 
for a cube of 200 h^ 1 Mpc size and a 256 3 mesh and for 
the same number of particles; each particle has a mass of 
4.0 x lO lo /i _1 M0. The transfer function was computed us- 



ing the CMBFAST code by Seljak and Zaldarriaga ( pL996| ). 
The rms mass density fluctuation parameter of this model 
is cr 8 = 0.87. 

We also calculated an additional model using the 
same code with 128 3 particles in a cube of size L = 
100 hr 1 Mpc, and the cosmological parameters fi m = 0.3, 
n A = 0.7, n b = 0.05, and h = 0.65; this COBE normalised 
model has the density fluctuation parameter eg = 0.78. 

We used the FoF algorithm to identify CDM halos. 
The only free parameter in the FoF method is the linking 
length, which is defined as the maximum separation be- 
tween particles which are still joined into groups. In our 
case the linking length was chosen as 0.23 in units of the 
mean particle separation, which approximately selects the 
matter density contrast Sn/n = 80. This density contrast 
was the one chosen by TUC to extract the group catalogue 
from LCRS, and it is typical of the values used in extract- 
ing groups from observational galaxy catalogues. While in 
the Einstein-de Sitter universe the overdensity within the 
virial radius of a cluster is 5n/n = 178, in a low density 
spatially flat universe with a cosmological constant it is 
Sn/n « 178ft-°- 6 (White et al. |1993[ ). The correspond- 
ing overdensity for our simulation is Sn/n ss 366. This 
means that the group sample determined using a low den- 
sity contrast (Sn/n — 80) contains groups which can be 
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in an uncertain dynamical state and do not have to be 
virialised. When particles outside the virialised core are 
included in the groups, the masses of simulated groups 
may be overestimated. For comparison we also calculated 
the mass function of publicly available numerical simula- 
tions of the Virgo consortium^, where Jenkins et al. 1998 
selected FoF-groups with a linking length of 0.2 of the 
mean inter-particle separation. 



3. Results 

3.1. Selection effects in the LCLG sample 

At first we discuss selection effects and show how we took 
them into account when calculating the mass function. 
There are two main selection effects in the LCRS: 1) ob- 
servations are performed in a fixed apparent magnitude 
interval which transforms to a distance dependent abso- 
lute magnitude interval, galaxies fainter or brighter than 
this interval are not included in the survey; 2) depending 
on the field 50 or 112 galaxies were measured for redshifts, 
the actual number of galaxies in the magnitude window 
could be larger, thus the sample is diluted. 

The first selection effect makes it impossible for galax- 
ies outside the window to enter the survey. Consequently 
groups consisting of faint galaxies occur only in the near- 
est region of the survey; with increasing distance fainter 
groups gradually disappear from the sample. This effect 
is seen in Figure |l| where the total luminosity of LCRS 
loose groups is shown as a function of distance. This se- 
lection effect can be statistically corrected using individual 
weights in calculation of the group mass function, follow- 



ing the procedure b y Mo ore, Frenk & White ( |l993|) and 
Girardi & Guiricin ( 2000 ). Namely, we weight each group 
by w = where T is the volume, which corresponds to 
the absolute magnitude limit of the third brightest galaxy 
of a group; if T > T (the full volume of a slice) , we take 
r = To. The data to calculate these limits have been 
tabulated by TUC separately for each of the 327 spec- 
troscopic fields of the LCRS. Figure || shows the distri- 
bution of the normalised weights wq (in units of the full 
volume of a slice) as a function of group mass. We see that 
there is a weak relation between cluster mass and the rel- 
ative weight: massive groups contain as a rule sufficiently 
bright galaxies and can be observed in the whole volume 
of the survey, thus for massive clusters the relative weight 
is closer to unity. 

Now we consider the second selection effect. To esti- 
mate the degree of dilution of the LCRS we compare the 
surface density of galaxies of the LCRS with the surface 
density of the Sloan Digital Sky Survey. Here the goal was 
to measure all galaxies in the magnitude window. The 
number of galaxies with measured redshifts in the Early 
Data Release of SDSS is 70 per square degree, while the 
LCRS has only 20 galaxies per square degree in the —6° 
slice, and 30 — 38 per square degree in the rest of slices. 
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Fig. 2. The weighting parameter wo as a function of the 
group mass 

The faint end limit in most slices of the LCRS is almost 
same as for the SDSS (—17.7), thus the difference in num- 
ber density is due to dilution of the LCRS. Most slices of 
the LCRS are diluted by a factor of about 2, and the slice 
—6° by a factor of 3.5. Dilution decreases not only the 
number of galaxies but also the number of clusters and 
groups, as the number of galaxies in a group may fall be- 
low 3 and the group will not be included in the catalogue. 
The number of groups in the —6° slice is actually lower 
than in other slices by a factor of about 2. 

However, uniform dilution does not influence strongly 
the mass function of groups. Firstly, virial masses of 
groups practically do not change. When determining virial 
masses, group members play the role of test particles, 
which move in a common gravitational potential. Using 
less test particles gives the same estimate for the poten- 
tial and for the virial mass as before, only with a larger 
variance. We checked this conclusion by diluting simulated 
group catalogues. 

Dilution may affect the group catalogue only by reduc- 
ing the group richness of smaller groups below the cata- 
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logue limit. These groups will drop out of the catalogue 
and the number density of smaller (less massive) groups 
will fall. We estimated this correction, using group cat- 
alogues obtained in our numerical simulations and ran- 
domly diluting the groups by 50%, as estimated above. 
We found that the conditional probability distribution 
P(log N\M) (N is the richness of a diluted group, M is the 
virial mass of the group, and we use decimal logarithms 
here) is close to a Gaussian with a mean 



log N{M) w log(M/M ) - 11.5 

and a rms error a that is also a function of the virial mass, 
a(M) « 0.3(14 - log(A//M Q )). 

The fraction a of the groups, which remain in the cata- 
logue after dilution, is given then by the error integral: 



a(M) = 0.5erfc 



logiVi 



log N(M) \ 
s/2a(M) J 



where Nl is the catalogue limit (3 and 5, in our case) and 
log N(M) and a(M) are given above. In order to restore 
the original differential mass function, we have to multiply 
it by the factor a -1 . This dilution correction is larger than 
1 only in the mass interval 10 12 -10 13 Mq, reaching 1.95 at 
its lowest end (for Nl = 3) and 2.67 for Nl = 5. As the 
differential mass function is small in this mass interval, 
this correction does not change appreciably the final inte- 
gral mass function. Also, the dilution correction does not 
change the high mass end of the mass function at all. 

3.2. Mass function of LCLGs 

The mass function (MF) of galaxy clusters/groups is de- 
fined as the number density of clusters above a given mass 
M , n(> M). To construct the group mass function from 
a group sample one needs accurate group masses and well 
defined volume and selection functions of the sample. We 
have used the masses estimated by TUC, who assumed 
that the groups were virialised and calculated virial masses 
of groups as: 



M*; ? > 



dtf os Rh 
G ' 



(1) 



where G is the gravitational constant, Rh is the harmonic 
radius of the group, oi os is the group line-of-sight velocity 
dispersion, and d = 6 in the case of a spherically symmet- 
ric velocity distribution of the group. 

We calculated the volumes of slices as follows: 



V 



cos S m AaA5 



4); 



(2) 



here 5 m is the mean declination of the slice, AS and Aa 
are the widths of the slice in declination and right ascen- 
sion (in radians), and r\ and ri are the lower and upper 
distance used in calculations. This formula is valid for spa- 
tially flat (k=0) cosmologies. 



Table 1. The numbers of groups in the LCLG catalogue. 
Columns give the mean declination of a slice, the number 
of all groups by TUC, the number of groups in the clean 
sample and the volume of the slice. 



Slice 5 


TUC 


Clean 


Volume [Mpc 3 /h 3 ] 


-3° 


288 


80 


927902 


-6° 


147 


37 


823452 


-12° 


276 


73 


874202 


-39° 


249 


71 


971817 


-42° 


257 


69 


988435 


-45° 


256 


63 


941531 


Total 


1473 


393 


5527339 


Mean 






921223 



The comoving cosmological distance r is a function of 
measured redshift and depends on the cosmological model 
as (e.g. Peebles 1993 ): 

dz' 

(3) 



H J E{z'Y 
where the function E(z) is given by 
E 2 (z) = n ro (l + zf + (1 - Q m - n A )(l + zf + fi A 



(4) 



The LCLG survey extends up to the redshift z ~ 0.15. 



Tucker et al. ( 2000 ) used the Einstein-de Sitter cosmology 
(assuming fl m = 1, fl\ — 0) to determine distances in the 
LCLG survey, while modern data prefer the parameters 
il m ~0.3 and fl/^—0.7. Thus for the present study we recal- 
culated comoving distances for all galaxies of the LCLG 
survey. For the particular redshift range 10000 < cz < 
45000 km s _1 the ratio between the comoving volumes in 
these two models is about 0.8. The numbers of groups for 
all slices, for the clean sample and the volume of each slice 
in the LCLG catalogue are shown in Table 1. Only groups 
with defined mass are included (for some groups the ve- 
locity dispersion, corrected for random errors, is zero and 
mass determination is impossible). 

The MF was constructed by sorting the masses in log- 
arithmic bins of width 0.2. Errors in group masses (an 
estimate of the rms error in virial mass has been found 
by TUC) were accounted for by convolving group masses 
with a log-normal distribution of the same mode and vari- 
ance. Figure || shows the mass function of LCLG for all six 
slices. To have a general picture we added slices together 
and found the average mass function, denoted by a thick 
solid line. We found that the average mass function and 
the mass function of the clean sample distributions are 
very similar, which indicates that the differences between 
the clean sample and the full sample are unimportant for 
the MF. We also tried excluding the slice S — —6° and 
found that it has only a weak influence on the total MF. 
In order to better show the differences between various 
slices we did not correct these mass functions for dilu- 
tion. The dilution-corrected final mass function is given 
in Table | and in Fig. g. 
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Fig. 3. The cumulative mass functions for all six slices 
in the LCLG. The thin curves from top to bottom are 
labelled as seen at the low mass end. The thick solid curve 
shows the average mass function for all the clusters. These 
mass functions are not corrected for dilution. 

Table 2. The cumulative mass function of LCLG 



log M [/i _1 M ] 


Density [h 3 Mpc~ 3 ] 


1.26E+12 


7.56E-04 


3.16E+12 


5.78E-04 


7.94E+12 


4.09E-04 


1.26E+13 


3.23E-04 


3.16E+13 


1.81E-04 


7.94E+13 


7.63E-05 


1.26E+14 


3.51E-05 


3.16E+14 


5.22E-06 


7.94E+14 


8.77E-07 



There is a clear difference in the MF between slices. 
The two northern slices at 5 = —6° and 5 = — 12° differ 
from others, having the lowest and the highest amplitudes 
at the faint end of the MF. We suppose that this discrep- 
ancy between mass functions are due to differences in di- 
lution, since the slice 5 = —6° was observed by the 50 fibre 
spectrograph only and the <5 = —12° slice by the 112 fibre 
spectrograph only. This discrepancy cannot be removed 
by a luminosity-weighting procedure. The remaining dif- 
ferences are probably due to cosmic variance. The slices 
of the LCRS survey intersect large voids and superclus- 
ters. In such slices, loose groups form a variety of struc- 
tures that may influence the observed mass function at its 
massive end, since the most massive groups (clusters) are 
associated with supercluster structures. 

Figure ^ shows clearly the flattening of mass functions 
at the low mass end. Such flattening is usually interpreted 
as the result of incompleteness of a sample (Girardi et al. 



199S), although it could also be real, describing a dimin- 
ishing number of groups of smaller masses. Thus, we have 
to study the completeness problem in more detail. 
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Fig. 4. The volume-density of groups as a function of dis- 
tance. Each line represents a different mass interval and 
is labeled by the decimal logarithm of its mean M/Mq. 



To estimate the mass completeness limit of the LCLG 
we can study the volume-density of the groups (within dif- 
ferent mass intervals) as a function of distance. In Figure |J 
lines denote mean values of the density in various mass 
ranges. Any rapid fall of the lines indicates incomplete- 
ness for that mass interval. We find from the Figure || that 
below masses of 10 135 — 10 14 Mq our sample is probably 
incomplete. 

Another way to estimate the mass completeness limit 
of the LCLG is to use the velocity dispersion distribu- 
tion function, VDF (Fadda et al. 



1996). The VDF is de- 



fined similarly to the MF, but for the velocity dispersion 
of groups. Assuming that the VDF can be described by 
a simple power law function (note that this assumption 
does n ot have a clear physical justification), Fadda et al. 
( 1996| ) set the completeness limit for their sample at the 
point where the power-law exponent of the VDF starts to 
change. For the LCLG this point is reached at the line-of- 
sight velocity dispersion er; os « 250 km/s. 

Figure || shows the relation between <ji os and mass M 
for the LCLG catalogue. Equation (|l|) defines the relation 
between these quantities via the harmonic radius of the 
group. The dotted line illustrates a power-law fit: 



hM 



10 15 M<r 







km/s, 



(5) 



where a = 819 and c = 0.35. The result is in quite good 
agreement with previous results by del Popolo & Gambera 
( |2000D and del Popolo et al. ( ^OOOQ , a = 842 km/s and 
c = 0.33, based on N-body simulations and X-ray obser- 
vations. A value 1/3 for the exponent b means that the 
system is virialised. The fact that the LCLG follow rather 
closely the virial relation shows that the possible contam- 
ination is small and the groups are physical. Setting the 
completeness limit to ai os = 250 km/s, based on the VDF, 
in equation (|5|), we get an estimate for the completeness 
limit in mass between 10 13 < M < lO 14 ft, _1 M , with the 
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Fig. 6. The left panel shows the velocity dispersion of dark matter particles in halos as a function of the halo mass. 
The one-to-one correspondence between the number of particles and the halo mass is adopted. The right panel gives 
the same relation, but unbound groups are removed. 
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for all groups and for groups with the number of members 
(N m > 5). This indicates that at least part of small groups 
with a large velocity dispersion (and thus a large mass) are 
probably pseudo-groups (studies by Einasto et al. 2002 in- 
dicate that such groups arc located mostly in low density 
regions). On the other hand, most of small groups have ve- 
locity dispersions in the same range as richer groups, and 
are probably real. In other words, exclusion of all small 
groups introduces an additional selection effect, which is 
difficult to quantify. 

4. Discussion 



1e+10 1e+11 1e+12 1e+13 1e+14 1e+15 4.1. Comparison with numerical simulations 



M[h"'M s 

Fig. 5. The velocity dispersion as a function of the mass 
for the LCLG catalogue. The line shows the analytical 
(least square) fit. 



mean value about M ~ 5 x 10 13 /i 1 Mq. This result agrees 
with the result of the first analysis. 



3.3. Reliability of groups 



Ramella et al. ( 1997 ) and Girardi & Giuricin ( p000[ ) noted 
that groups with a few (less than 5 members) may consti- 
tute not real, but pseudo-groups. On average, in pseudo- 
groups the velocity dispersion is considerably larger than 
in real groups (Ramella et al. 1997). Indeed in the LCLG 



sample, groups with less than 5 galaxies have a very large 
scatter of velocity dispersions, and some velocity disper- 
sions are extremely large. Analysis of a subsample, where 
poor groups with only 3 or 4 members arc excluded, shows 
that in the w — M plane (as in Figure ||) the scatter for this 
subsample is somewhat smaller than for the whole sam- 
ple. However, exclusion of all poor groups decreases the 
mass function considerably in the whole mass range (see 
Figure [?], where the mass function is given for two cases, 



Now we compare our mass function with mass functions 
found for groups in numerical simulations. Figure^ shows 
the velocity dispersion of dark matter particles in groups 
as a function of the halo (group) mass in our numerical 
simulations. We see that in small groups with mass less 
than ~ 10 13 /i _1 Mq (which contain less than ~ 250 parti- 
cles) the scatter of the velocity dispersion is very large. 
This phenomenon can be explained if we assume that 
small groups are not yet virialised and are contaminated 
by non- virialised particles, or alternatively are composed 
of several small subgroups moving fast relative to each 
other. Another explanation for this effect is that it is a 
result of the chosen 8n/n — the FoF algorithm collects 
particles from outskirts of small clusters although they 
do not belong dynamically to the group. The fraction of 
interlopers in simulated loose groups can become rather 
high, about 20-30%, as found by Diaferio et al. ( |1999| ). 

To test the second assumption we built groups, using a 
lower linking length b = 0.17 (in units of the mean particle 
separation). This did not change considerably the scatter 
of the velocity dispersion among the groups with less than 
250 particles. Thus the first explanation is more likely. 

To eliminate the effect of unbound groups in simula- 
tions we calculated both the potential and kinetic energy 
for all groups and excluded unbound particles, using the 
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condition \E pot \ < 0.8Eki n (E pot is the potential energy 
and Ekin the kinetic energy of a group) for an unbound 
particle. We split groups with a high kinetic energy as 
compared to the potential energy to separate interact- 
ing or merging groups. To this aim we employed a 6D- 
group finder with a linking measure Ar 2 /r 2 + Ai^/cr 2 , 
where r and a v are the group effective (half-mass) radius 
and the velocity dispersion. This procedure affected only 
about 10% of the halos, with many of the halos having a 
small number of DM particles. We show the results of our 
clean-up procedure in the right panel of Figure [| which 
is similar to the left panel, only non-virialised particles 
and bypassing groups are excluded. We see that our pro- 
cedure to eliminate non-virialised particles works well — 
the low-mass groups follow the same trend as observed for 
the high-mass groups and clusters. 

We have calculated the mass functions for the sample 
of all groups and for the sample of virialised groups in 
our simulations. These functions are very similar, only for 
groups of very small mass (M < 10 12 /i _1 Mq) the MF 
of virialised groups differs from that of all groups by less 
than 1%. Thus the contribution of unbound groups to the 
total mass function is practically negligible. 

Figure shows the mass function of dark matter halos 
in our simulation together with the MF of the LCLG. The 
Poisson error bars (square root of counts per bin multi- 
plied with the mean weight) for the full LCLG mass func- 
tion are shown by the dash-dotted lines. In order not to 
overcrowd the figure, we have not shown the error bars for 
the N >— 5 sample; they are very close to those of the 
full sample. For small masses (M < 5 • 10 13 Mq) the error 
bars are smaller than the size of the plotting symbols. For 
comparison, Figure shows also the mass function of the 
simul ations by the Virgo Cons ortium , ACDM-gif (Jenkins 
et al. 1998 , Kauffmann et al. 1999a ). In the Virgo simu- 
lation the side of the box was 141.3 Mpc, the matter 
density O m = 0.3, the cosmological term J7a = 0.7, the 
particle mass 1.4 x 10 10 Mq//i, and the number of the 
particles 256 3 . Figure shows that our 128 3 simulation 
with the amplitude parameter er 8 = 0.78 gives the best fit 
among the models we used to simulate the LCLG subsam- 
ple at the massive end (2 • 10 14 M Q - 1O 15 M ) of the MF. 
For smaller masses, down to 4 • 10 13 Mq the model with a 
higher cr 8 = 0.87 gives a better fit. The cluster MFs in the 
Virgo Consortium simulation and in our 256 3 simulation 
are quite similar. For yet smaller masses, where we sus- 
pect that the LCLG sample is incomplete, the LCLG mass 
function flattens out, in contrast to the rising simulated 
mass functions. 

Most of earlier results give er 8 ~ 0.9—1 for the 
fl m = 0.3 scenario (Eke et al. 1996, Viana and Liddle 1996 



etc.). This value agrees with the recent result by Pierpaoli 
et al. (2001), <r 8 == 1-02. Lower values are obtained based 



on X-ray cluster samples. Reiprich & Bohringer (2002) 
derived for ft m = 0.3, cr 8 - 0.68. From the ROSAT Deep 
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Fig. 7. Comparison of the dilution-corrected observed 
mass functions and the theoretical mass functions. The 
observed mass functions are shown with symbols: filled 
circles (the LCLG whole sample), open circle s (the LCLG 
N m > 5) and stars (Girardi & Giuricin (|2000|) . The dashed 
line shows our 256 3 simulation MF (erg = 0.87), the solid 
line shows our 128 3 simulation (ct 8 = 0.78) MF, and the 
dotted line shows the Virgo Consortium simulation, MF. 
For clarity Poissonian errors are shown only for one case 
(the full mass function). 



et al. ( P002[) used the SDSS data and got <r 8 = 0.72 for 
&m = 0.3. In principle, <r 8 can be determined on the basis 
of the CMB data only; future experiments as MAP and 
PLANCK will help to pinpoint this parameter. 

4.2. Comparison with earlier work 

The mass function determined by Girardi & Giuricin 
( 2000| ) is shown in Figure |. This group catalogue was 
constructed using a subsample of galaxies of the Lyon- 
Meudon Extragalactic Database. Groups were identified 
using two different methods: the FoF algorithm (P-groups, 
453 groups) and the hierarchical method (H-groups, 498 
groups). These group samples extend until cz < 5500 
km s _1 . 

Here we used the P-data from figure 3 in Girardi & 



Cluster Survey Borgani et al. (2001) derived Cl m — 0.35 
and erg ~ 0.66. Recently Lahav et al. (2002) obtained 
cr 8 = 0.73 using the 2dFGRS+CMB data. Finally, Bahcall 



Giuricin (2000). In that case they chose the subsample of 
the nearby groups with cz < 2000 km s _1 . , which they as- 
sumed to be a complete sample and a good representation 
of the total population. However, this sample covers only 
the local supercluster, while the LCLG catalogue extends 
up to distances ~ 450 h" 1 Mpc. This volume contains al- 
ready a variety of structures. 

We see that in the high mass range (M > 10 14 /i _1 Mq) 
the MFs are rather similar. Around masses M w 5 • 10 13 
the mass functions of the simulations and the LCLG start 
to differ from each other. Also, at t hese lo wer masses the 
mass function by Girardi & Giuricin ( 2000| ) is considerably 
higher than the LCLG mass function. 

This is an intriguing result that certainly needs fur- 
ther analysis. We see that although loose groups promise 
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to provide better statistics, there are a number of un- 
solved problems that do not permit a realisation of this 
promise at the moment. There are two problems in the 
LCRS group sample: reliable determination of masses of 
small groups, and selection effects. Virial mass estimates 
can be in error; Girardi and Giuricin (200C) try to correct 
for this. Their correction is rather complicated and could 
induce its own errors. 

The simplest explanation of the difference between the 
mass functions found in the present paper and that of 



Girardi and Giuricin (200C) is a mass-dependent selection 
effect. The LCRS is a strongly diluted sample of galaxies, 
and it is biased against low surface brightness galaxies. 
However, as shown above (Sec. 3.2), dilution practically 
does not affect the integral mass function. 

The second effect, the absence of low-brightness galax- 
ies is probably much less present in the Lyon-Meudon 
galaxy catalogue, from which Garcia ( 1998 ) extracted the 
catalogue of local grou ps th at served as a source for the 
Girardi and Giuricin ( [2000 ) mass function. If low sur- 
face brightness galaxies were more prevalent in low mass 
groups - which is plausible - the LCLG mass function 
would flatten out as seen in Fig. |^. 

In a recent paper Martinez et al. ( |2002 ) found a similar 
flattening of the mass function of groups, extracted from 
the 2dF catalogue. This is also a deep sample, similar to 
the LCLG catalogue. 

We should also look for possible problems with our 
simulated mass functions. 

One problem is that we are comparing dark matter 
halos with galaxy-populated halos. It is probable that ha- 
los of smaller masses (10 12 — 10 13 ) can not easily host 
three or more galaxies. This would explain the deficiency 
of small-mass galaxy groups. Thus, better simulations of 
group catalogues will have to include galaxy formation. 

There could also be a problem with mass determina- 
tion. We devoted much attention to accurately evaluate 
the mass estimates in simulations; in particular, we suc- 
ceeded in filtering out unbound groups. Still, the mass- 
velocity dispersion relation of virialised groups widens to- 
ward the small mass end, showing that virial masses of 
small groups may have large intrinsic errors. This widen- 
ing could be caused, first, by the predominance of radial 
orbits in small groups, but this should be a rather small 
effect. Another possibility is that small groups might have 
a higher intrinsic velocity spread than large groups due to 
different formation histories, or in other words, due to a 
mixture of young and old groups. This is a question that 
could be only answered by a detailed dynamical study of 
observed groups. 

Another numerical effect could be an enhanced group 
production in simulations, caused by gravitational two- 
body collisions, which were discussed by Suisalu & Saar 



1995 and Splinter et al. 1998). More recent analysis con 



5. Conclusion 

In this work we have presented the mass function of the 
Las C ampanas Loose Groups selected by Tucker et al. 
( 2000| ) . We have also studied the mass- velocity dispersion 
relation that shows the expected scaling relation for viri- 
alised systems. We conclude that the loose galaxy groups 
are basically physical systems. The mass function was 
compared with results of dark matter simulations. 
Our main conclusions are: 

1) The completeness interval of the LCLG sample is 
found to be around 5 • 10 13 M Q //i < M < 8 x 10 14 Mq//i. 

2) The high end of the mass function of the LCLG 
sample lies between the simulated mass functions corre- 
sponding to erg = 0.78 and as = 0.87. 

3) At small masses M < 5 • 10 13 Mq//i the mass func- 
tion of the LCRS sample of loose groups flattens out. 
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